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A study of the tetragonal to orthorhombic phase transition of GdFeAsO is presented. Planes of 
the reciprocal space were reconstructed form single- crystal synchrotron X-ray diffraction data. By 
cooling below the structural transition temperature splitting of the Bragg reflections was observed 
corresponding to four different twin domain orientations. A model was developed to quantify the dis- 
tortion of the lattice from the position of the splitted reflections relative to each other. Constrained 
2D-Cauchy fits of several splitted reflections provided positions of the reflections. The influence of 
the structural distortion was detectable already above the structural transition temperature hinting 
at fluctuations in the tetragonal phase. 
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I. INTRODUCTION 

With the discovery of high-temperature supercon- 
ductivity in layered iron containing compounds (iron 
pnictides , a new system of phases was added to the 
study of unconventional superconductivity. Since the un- 
derstanding of the phenomenon of high-T c superconduc- 
tivity is still limited, the question of the differences and 
similarities towards other high-T c compounds such as the 
cuprates arises*. 

The cuprates as well as many iron pnictides become 
superconducting by doping of a non-superconducting an- 
tiferromagnetic "parent" compound. Whereas the su- 
perconductivity in the cuprates is believed to originate 
from the Mott-insulator properties of the "parent" com- 
pound, the iron pnictides do not exhibit this electronic 
correlation. However, a common feature is the antiferro- 
magnetic ordering, which has to be suppressed to achieve 
superconductivity and is thus considered significant. 

Iron pnictides, such as REFeAsO, (RE = rare-earth 
metal) and ^4Fe2As2 (A — Ca, Sr, Ba, and Eu), show 
antiferromagnetic transitions accompanied by structural 
distortions from tetragonal to orthorhombic symmetry 
upon cooling. While the structrual transition and the 
magnetic ordering occur simultaneously in AF^Asu^, 
the structural transition in REFeAsO compounds is ob- 
served at slightly higher temperatures than the magnetic 
ordering^— . 

The i?EFeAsO room-temperature phases of ZrCuSiAs- 
type structure (space-group type P4/nmm) undergo 
a translationengleiche transition of index 2 to the low- 
temperature structure with space-group type C mmeZ- 
(Fig. [TJ for structural relations see also Ref. [l2T ). A 
similar transition occurs in the AFe2As2 compounds 
(J 4/ 'm mm to F mm m). Tanatar et al*£ and Blombcrg 
et alr^- studied the formation of lamellar transformation 
twins by single-crystal X-ray diffraction of AFe2As2 and 
the influence of the domain structure on the anisotropy 




FIG. 1. (Color online) Low-temperature crystal structure of 
GdFeAsO at 50 K (space group Cm me, No. 67). The dis- 
placement ellipsoids represent 99.9 % probability. Solid lines 
between iron atoms indicate shortest Fe-Fe distance, dashed 
lines are slightly longer, ha is the height of gadolinium ions 
above the oxygen layer, /ia s is the height of the arsenic ions 
above the iron layer, and u) c is the As-Fe-As angle with the 
angle bisector parallel to the c-direction. 



of the resistivity in directions of the basal plane. Due to 
the analogous mechanism of the structural transition for 
REFeAsO, the question arises, if the twin formation is 
similar to /!Fe2As2. 

Here, we present the first single-crystal X-ray diffrac- 
tion study of the structural evolution and the twin for- 
mation of REFeAsO compounds with GdFeAsO as rep- 
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resentative. Doped GdFeAsO has one of the highest T c 
(56 K) in this class of superconductors^, which has been 
associated with the vanishing distortion of the FeAs4 
tctrahedra 16 . This concept has been the starting point of 
an ongoing investigation on the structural changes asso- 
ciated with the transition to the superconducting phase 
and its physical implications. 



II. EXPERIMENTAL 



Preparation and characterization of the crystal used for 
single-crystal synchrotron X-ray diffraction has been de- 
scribed elsewhere^. Diffraction data for this study were 
recorded with a 165 mm MAR-CCD detector mounted 
on a Huber four-circle diffractometer at beamline D3 at 
DESY. The measurements were performed in the tem- 
perature range of 50-300 K using an open flow Oxford 
Diffraction Hclijet cryostat. </?-scans were recorded with 
an increment of 1°, a detector distance of 54.4 mm, and a 
wavelength of 0.49741 A. The detector distance and the 
wavelength were refined by comparing the experimental 
lattice parameters of a standard corundum single-crystal 
to literature data. The raw frames from the MAR-CCD 
detector were converted using the APEXII suited. Inte- 
gration and corrections for oblique incidence and polar- 
ization were performed within SAINT 4^. For data re- 
duction and absorption correction SADABS 2 " was used. 
The structures were solved and refined with SHELXS and 
SHELXL2i. Details concerning the structure analysis of 
GdFeAsO at different temperatures can be obtained from 
the Fachinformationszentrum Karlsruhe 2 ^. 



III. RESULTS AND DISCUSSION 

When cooling below the temperature of the structural 
transition, T s sa 115 K, splitting of reflections at high 
diffraction angles was observed in the synchrotron exper- 
iment. To assess the metric distortion during the struc- 
tural transition, a program was written in MATLAB 23 , 
that allows reconstruction of planes of the reciprocal 
space with high lateral resolution (Fig. ^ . 

The metric distortion of the structure is the only de- 
gree of freedom generated by the structural transition. 
A model has been developed, which allows to determine 
the metrical distortion from the observed splitting of re- 
flections. For convinient comparison of the tetragonal 
(indicated by subscript t) and orthorhombic (subscipt 
o) structures, the C-centered orthorhombic cell is trans- 
formed into a monoclinic (subscript m) primitive cell (up- 
per right panel in Fig. [3]) by 

a m = ^ (a Q + b Q ) and b m = ^ (-a c + b Q ) 



with |a n 



and laJ > |b Q 




FIG. 2. Reconstructed h k plane in monoclinic indexing of 
GdFeAsO at 50 K. Inset: Magnification of the 6 4 reflections. 



The twin domains join in common planes that are in- 
dexed (1 1 0) and (1 1 0) in the C-centered cell, i.e. (1 0) 
and (0 1 0) in the primitive cell. Thus four different ori- 
entations of the domains in the low-temperature phase 
result. In the diffraction pattern fourfold splitting of re- 
flections in general positions can be observed, especially 
at higher diffraction angles (Inset in Fig. [5]). The position 
vectors pi, P2, P3, and P4 of the reflections belonging to 
the four domains can be derived from their common cen- 
ter p c and the distortion angle 7* (Fig. [3]) by 



Pi 

P3 



1 -sin( 7 * - 90°) 
cos( 7 * - 90°) 

cos( 7 * - 90°) 
- sin( 7 * - 90°) 1 



Pc, P2 



Pc, P4 



1 sin( 7 * - 90°) N 
cos( 7 * - 90°),' Pc 

/cos( 7 * - 90°) 0^ 
^sin( 7 * - 90°) \ t 



The distortion of the structure is then defined by the 
angle 7 > 90° (i.e. 7* < 90° in reciprocal space). 



Blomberg et al— detwinned BaFe2As2 crystals by ten- 
sile stress and showed that the relative positions of the 
reflections is dependent on the domain fractions. If, for 
example, the domain fraction of D4 in Fig. [3] increases, 
the orientation of the two sets of domains changes to- 
wards each other. In the geometrical model described 
here, the correlation of the reflection positions to their 
intensity is achieved by linear scaling of the transforma- 
tions shown above. This approximation is acceptable due 
to the small differences of the intensities and the local fit 
of the reflections. In summary, the positions of the re- 
flections are constrained to their common center p c of all 
the four reflections, the monoclinic distortion angle 7*, 
and the intensities of the reflections. This enables the 
parameter refinement of overlapping reflections. 

To quantify the distortion, an algorithm was written in 
MATLAB to fit the splitted reflections by a linear com- 
bination of four elliptical 2D-Lorentz (Cauchy) functions 
(concerning multivariate t-distributions see Ref. 24, 2D- 
Gaussian fit of neutron data see Refs. [25| and 26). A 
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FIG. 3. (Color online) Scheme of the tetragonal (left) to orthorhombic (right) distortion in analogy to Ref. [i~3l Upper panels: 
Schematic reflection patterns for a* fe*-planes with I 7^ 0. Lower panels: Real space domains. Black lines represent the primitive 
unit cells, green points are the positions of iron atoms and green lines show shortest Fe-Fe distance. Dl to D4 denote the domain 
orientations corresponding to pi to p4 in the upper right panel. 




FIG. 4. (Color online) Contour plot of the 5 50 reflection 
in hkO plane. Inset: ID-section of the 2D-distribution as 
indicated by the black line in the contour plot. 

Lorentzian type profile function was chosen empirically 
by investigating ID-sections of the 2D-distribution (see 
Fig. [4]). The fitting parameters were the center of the 
splittccl reflections p c , the distortion angle (7*), the el- 
lipticity of the reflections, the width of the 2D-Lorentz 
distributions, and the individual intensities of the reflec- 
tions. 

For the structural investigation 20 reflections of each 
measurement, exhibiting reflection splitting, were fitted 



to obtain an average distortion angle 7* and a corre- 
sponding standard deviation. The reflection data were 
integrated with a tetragonal cell since splitting was not 
wide enough to perform multidomain integration. Using 
the tetragonal lattice parameters, which correspond to 
Pc and the obtained distortion angle, the orthorhombic 
lattice parameters were obtained by 

d = 2 d t sin ^— ^ and b Q — 2 a t cos ^— ^ . 

All reflections could be fitted with this model. No signs 
for a deviation from the low-temperature orthorhom- 
bic structure towards true monoclinic symmetry was ob- 
served. Furthermore, no superstructure reflections could 
be identified in the reconstructed images of the recipro- 
cal space. This confirms the tetragonal to orthorhombic 
transition established by powder diffraction^. 

The order parameter of the tetragonal to orthorhombic 
transition is proportional to the spontaneous deformation 
e and can be fitted by a power law (upper panel in Fig. [5j, 

a~b fj-90°\ , / TV 

e = = tan = A 1 

a+b \ 2 J \ Tj 

with the critical exponent j3 = 0.072(1) and a transition 
temperature of T s = 111(1) K. 

The critical exponent is very small but comparable to 
the order parameters found for SrFe2As2 and EuFe2As2 
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FIG. 5. (Color online) Temperature dependence of selected 
structural parameters of GdFeAsO single crystal. Error bars 
represent ±<r. Upper panel: Spontaneous deformation fit- 
ted by a power law. Middle panel: Deviation of the lat- 
tice parameters relative to 300 K. (a' t = a t /a t (300K) — 1, 
a' a = Oo/ [V2a t (300K)] - l,b' = b a / [y/2 a t (300 K)] - 1, and 
c' = c/c(300K) — 1). Lower Panel: Cell volume of the or- 
thorhombic phase and doubled volume of the tetragonal phase 
for comparison. Lines are guides for the eye in the middle and 
lower panel. 



(0.098(1) and 0.112(1), respectively), where structural 
and magnetic transition occur at the same temperature. 
Tegel et al. explained the small critical exponent by an 
ordering following the two-dimensional Ising model of the 
iron subsystem, which should yield a critical exponent of 
1/8. In the case of GdFeAsO the structural transition 
occurs at temperatures above the magnetic transition, 
which could be the reason for the even lower critical ex- 
ponent determined, indicating that the structural transi- 
tion of GdFeAsO is close to first order (concerning mul- 
tiferroic transformation and the connection to first order 
phase transition see also Refs. andl28r). 

The structural transition temperature T s = 111(1) K 
obtained from the fit is lower than determined by Luo 
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FIG. 6. (Color online) Temperature dependence of the As- 
Fe-As angles oj with its bisector parallel to the direction in 
subscript. Note the relation cosaj at = — ~ (1 + cosoj c ). Lines 
are guides for the eye, error bars represent ±<r. 



et al£ with powder X-ray diffraction and resistivity mea- 
surements. This deviation might be caused by a system- 
atic temperature error from the open flow helium cryostat 
used for this single- crystal X-ray diffraction study. 

The lattice parameters and the volume of the tetrag- 
onal phase (middle and lower panel of Fig. [5]) decrease 
in the course of cooling. However, there are hints that 
there is an anomaly just above the structural transition 
(e.g. different slope for c below 140 K), as it was also 
observed for the thermal expansion of GdFeAsO by Klin- 
geler et alJ^- and assigned to a competition of differently 
ordered phases. 

Interestingly, the height of the gadolinium atom above 
the oxygen layer (/led in Fig.[T|) remains almost constant 
over the whole temperature range, while the arsenic atom 
height over the iron layer (/lAs) diminishes with decreas- 
ing lattice parameter c. The angle w c of the coordination 
polyhedron around Fe (Fig. [5]) with its bisector parallel 
to the c-direction first increases upon cooling and then 
becomes almost constant already above the structural 
transition and does not change much upon further cool- 
ing. This indicates structural fluctuations by softening 
of the lattice or vibrational modes above the structural 
distortion which is also mirrored in other properties (e.g. 
thermal expansion^ and X-ray reflection broadening^). 
The angles with the angle bisector parallel a Q and b Q do 
not exhibit this anomaly and merely differentiate corre- 
sponding to the tetragonal to orthorhombic transition. 

Although, the overall structural changes upon cooling 
are subtle, they indicate a response of the structure al- 
ready at temperatures above the onset of the metric dis- 
tortion. It is clear that higher resolution of the structural 
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investigations is needed to further study the interplay of 
structure and physical properties of both the "parent" 
and superconducting phases. 

IV. CONCLUSION 

A geometrical model was developed for single- crystal 
synchrotron X-ray diffraction measurements to quantify 
the metric distortion and twin formation of GdFeAsO in- 
duced by the tetragonal to orthorhombic structural tran- 
sition. 

Subtle changes in the structure are detectable. Fit- 
ting the spontaneous deformation by a power law lead 
to a small critical exponent of 0.072(1), which is lower 
than expected for a two-dimensional Ising model, point- 
ing towards a transition close to first order. This may be 
caused by the separation of the magnetic and the struc- 
tural transition in temperature. The lattice parameters 
decrease with decreasing temperature, and the a and b 
axes differentiate at temperatures below the structural 
transition. The deformation is abrupt and preceeded by 
a stagnation of the decrease of the cell parameters. This 



indicates long-range fluctuations of the lattice or soft- 
ening of a vibrational mode at temperatures above the 
structural transition, which were also observed by other 
methods (e.g. Ref. l30j). If this feature is linked to the 
superconductivity is not yet clear. 

Although the presented structural information from 
single-crystal X-ray diffraction show the subtle structural 
changes around the tetragonal to orthorhombic phase 
transition, data with higher resolution is needed to un- 
reavel more details of the interplay of the structure and 
properties of the "parent" compounds of iron-based su- 
perconductivity. 
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